First-forbidden mirror /3— decays in ^4 = 17 mass region 



N. Michelf, J. Okolowiczt§, F. Nowacki|, and M. Ploszajczakf© 
t Grand Accelerateur National d'lons Lourds (GANIL), CEA/DSM - CNRS/IN2P3, BP 5021, 

F- 14-076 Caen Cedex 05, France 
§ Institute of Nuclear Physics, Radzikowskiego 152, PL - 31342 Krakow, Poland 
:(: Laboratoire de Physique Theorique Strasbourg (EP 106), 3-5 rue de VUniversite, F-67084 

Strasbourg Gedex, France 
©Institute for Nuclear Theory, University of Washington, Box 351550, Seattle, WA 98195, USA 



Abstract 



The first-forbidden /J— decay of "'^'^Ne into the 'halo' state J'^ = '\-/2\ of ^''F presents 
one of the largest measured asymmetries for mirror /?— decay feeding bound final 
states. This asymmetry is studied in the framework of the Shell Model Embedded 
in the Continuum (SMEC). The spatial extent of single particle orbits is constrained 
by the proton capture cross-section ^^0(p, 7)^'^F calculated in SMEC. This allows 
to estimate mirror symmetry breaking in ^''F/^^0 and ^''Ne/^''N nuclei. 

21.60. Cs, 23.40.-S, 23.40.Hc, 25.40.Lw 
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I. INTRODUCTION 



A realistic account of the low-lying states properties in exotic nuclei requires taking into 
account the coupling between discrete and continuum states which is responsible for unusual 
spatial features of these nuclei. This aspect is particularly important in the studies near drip 
line where the amount of experimental information is strongly reduced and one has to use 
both structure and reaction data to understand basic properties of these nuclei. Within 
the Shell Model Embedded in the Continuum (SMEC) approach one may obtain the 
unified description of the divergent characteristics such as the spectra (energies of states, 
transition probabilities, proton/neutron emission widths, /5— decays, etc.) and the reactions 
involving one-nucleon in the continuum (proton/neutron capture processes. Coulomb disso- 
ciation reactions, elastic/inelastic proton/neutron reactions, etc.). This provides a stringent 
test of the effective interactions in the SMEC calculations and permits to asses the mutual 
complementarity of reaction and structure data for the understanding of these nuclei. In 
this context, it can be interesting to compare the first-forbidden /3— decay transition of ^^Ne 
in the ground state (g.s.) J'^ = l/2f to the weakly bound, first excited state J'^ = 1/2^ in 
^^F, with the corresponding mirror decay of ^^N into a well bound excited state of ^''O. In 
these decays an abnormal asymmetry of mirror decay rates has been observed by Borge et 
al. and later confirmed by Ozawa et al. [Q. Borge et al. explained this effect by the 
large asymmetry of radial sizes of Si/2 s.p. orbits involved in bound states of ^''F/^''0 and 
^^Ne/^^N. Different explanation has been provided by Millener |^ who attributes the bulk 
of the asymmetry to charge dependent effects which lead to different lsi/2 occupancy for the 
initial states , i.e. to different amplitudes of 7i{0pl^2^sl^2)^i^Pi/2) ^i^Pi/2^^i/2)''^i^Pi/2) 
components in the g.s. of ^^Ne and ^^N, respectively. 

In the framework of SMEC, we shall discuss the constraints from proton capture cross- 
section data on the radial wave function involved in the description of the first-forbidden 

decay into bound final states. 

In the SMEC formalism, which derives from the continuum shell model formalism [^0, 
the subspaces of (quasi-) bound (the Q subspace) and scattering (the P subspace) states are 
separated using the projection operator technique. P subspace contains asymptotic channels, 
which are made of (A^ — l)-particle localized states and one nucleon in the scattering state, 
whereas Q subspace contains many-body locahzed states which are build up by the bound 
single-particle (s.p.) wave functions and by the s.p. resonance wave functions. The wave 
functions in Q and P are then properly renormalized in order to ensure the orthogonality 
of wave functions in both subspaces. 

In the first step, we calculate the (quasi-) bound many-body states in Q subspace. For 
that we solve the multiconfigurational Shell Model (SM) problem : Hqq^i = Ei^i , where 
Hqq is given by the SM effective interaction which is responsible for the internal mixing of 
many-body configurations. The quality of the SMEC description depends crucially on the 
realistic account of configuration mixing for the coexisting low-lying structures and hence on 
the quality of the SM effective interactions and the SM space considered. In the description 
of A = 17 nuclei, it is important to take into account the dynamics of ^^O core and to 
include 2p-2h and 4p-4h excitation from p— shell to sd— shell. Zuker-Buck-McGrory (ZBM) 
effective SM interaction in the basis of Opi/2, lsi/2 and 0^5/2 orbitals allows to take 

into account this dynamics. The valence space (Opi/2, lsi/2, 0^5/2) has the advantage to be 
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practically non spurious and most of states at the p — sd interface are well described through 
the configuration mixing of these three orbitals. Results of this paper have been obtained 
with the ZBM interaction. 



II. DETAILS OF THE MODEL 



To generate the radial s.p. wave functions in Q subspace and the scattering wave functions 
in P subspace, as a first guess, we use the average potential of Woods-Saxon (WS) type with 
the spin-orbit : Vso^l-i'2l-s)r~^ df (r) / dr , and Coulomb parts included. = 2fm^ is the pion 
Compton wavelength and /(r) is the spherically symmetrical WS formfactor. The Coulomb 
potential Vc is calculated for the uniformly charged sphere of radius Rq (see Table 1). This 
'first guess' potential U{r), is then modified by the residual interaction. We shall return to 
this problem below. 

For the continuum part, we solve the coupled channel equations : 

{E - Hpp)i^i~^ ^ j:iE - H^c')d^-'^ = , (1) 

c 

where index c denotes different channels and Hpp = PHP. The superscript (+) means 
that boundary conditions for incoming wave in the channel c and outgoing scattering waves 
in all channels are used. The channel states are defined by coupling of one nucleon in the 
scattering continuum to the many-body SM state in (A^ — l)-nucleus. For the coupling 
between bound and scattering states around ^^O, we use the density dependent interaction 



(DDSMl) pO| , p!l| . This interaction provides external mixing of SM configurations via the 
virtual excitations of particles to the continuum states. The channel - channel coupling 
potential : H^^r = (T + U)6^' + , , contains the kinetic-energy operator T and the 
channel-channel coupling v^^, generated by the residual interaction. At a first step, the 
potential for channel c consists of the initial WS potential U{r), and of the diagonal part of 
coupling potential v;!^ which depends on both the s.p. orbit (pij and the considered many- 
body state J. Hence, the initial potential is modified by the couphng potential and in the 
next step the s.p. wave functions (pij defining Q subspace are generated by the modified 
potential, what in turn modifies the diagonal part of the residual force, etc. In other words, 
the procedure of solving coupled channel equations (|I]) is accompanied by the self-consistent 
iterative procedure which for each total J yields the J— dependent self-consistent potential 
+ v^^^'^'' (r) , and consistent with it the new renormalized formfactor of 
the coupling force. U^^^\r) differs significantly from the initial potential, especially in the 
interior of the potential For weakly bound many-body states, strong modification of 

the surface features of the initial potential U{r) has been found as well P,p!0|. Parameters 



of the first guess potential U{r) are such that U^'^'^^r) reproduces energies of experimental 
s.p. states, whenever their identification is possible. For certain J, the s.p. wave functions 
are not modified by the above iterative procedure. For example in J'^ = 1/2+ states of ^''F 
and ^''O, the couplings to the g.s. 0+ of ^^O modify only lsi/2 s.p. wave function. In this 
case, the radial wave functions of all other s.p. states are generated by an auxiliary average 
potential [/(""^■)(r) of the WS type with spin-orbit and Coulomb parts included. We shall 
return to this problem below. 
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The third system of equations in SMEC consists of the inhomogeneous coupled channel 
equations: 

{E^+^ -Hpp)ujt^ = HpQ<l>, = Wi . (2) 

The source term Wi depends on the structure of iV - particle SM wave function $j. The 
radial formfactor of the source depends on s.p. wave functions of U{r) (U^"''^^\r)). Solutions 
of the eqs. @ describe the decay of quasi-bound state $i in the continuum. Reduced 
matrix elements of the source term involve products of two annihilation operators and one 
creation operator : 'R--'^s{l)i3 ~ {(^^i3{(^'y(^s)^y°' ■ In the SMEC calculations for ^^F and ^''O, the 
matrix elements of the source term depend sensitively on the percentage of the shell closure 
in ^^O, i.e., on the amount of correlations both in the g.s. of ^^O and in the considered states 
of ^^F or ^^O [O. The total wave function is expressed by three functions: $j , 



where : 



n = CE + + I Hqp \ ^e) (3) 



H'^^{E)=Hqq + HqpG^^^Hpq , (4) 



is the energy dependent effective SM Hamiltonian in Q subspace which contains couplings 
to the continuum. Operator Hqq{E) is hermitian for energies below the particle emission 

threshold and non-hermitian for energies higher than the threshold. The eigenvalues Ei—^iTi 
are complex for decaying states and depend on the energy E of the particle in the continuum. 
The energy and the width of resonance states are determined by the condition: Ei{E) = 
E. The eigenstates corresponding to these eigenvalues can be obtained by the orthogonal 
but in general non-unitary transformation [jl],|^,|^J^ . Inserting them in one obtains 
symmetrically the new continuum many-body wave function modified by the discrete states 



i E- Ei + {i/2}Ti 

and the new discrete state wave function modified by the coupling to the continuum states: 

n^ = ^^ + J: dE Ce' ^(+) _ ^, {Ce' I H I $.) . (6) 

These SMEC wave functions will be used in this paper to calculate spectra of ^''F, ^'^O, 
the first-forbidden /3— decays : ^^Ne(/5"'')^'^F, ^'^N(/3~)^''0, and the radiative proton capture 
reaction -'^^0(p, 7)"'^'^F. 

Figs. 1 and 2 show SMEC energies and widths for positive parity (l.h.s. of the plot) 
and negative parity (r.h.s. of the plot) states of ^^F and ^^O, respectively. Large breaking 
of mirror symmetry, which can be seen by comparing spectra in Figs. 1 and 2, is due to 
the combined effect of the low separation energy and the Coulomb field. The continuum 
coupling, which due to different positions of the lowest particle emission thresholds acts 



4 



differently in ^^F and ^^O, cannot fully account for the observed mirror symmetry breaking. 
The simplest way of correcting this deficiency is to adjust the spacing of 6/5/2 and Si/2 s.p. 
orbitals in ZBM interaction in such a way that the experimental energy splitting between 
the g.s. 5/25*" and the first excited state 1/2^ in ^'^O and ^'^F is reproduced by the SMEC 
calculation. In this way, the s.p. energy of ^5/2 in ZBM interaction becomes : Ed^^^ = 3.95 and 
3.5 in ^''F and ^''O, respectively. The energy of Si/2 orbital remains £si/2 = 3.3 in both cases. 
In the following, these two hybrid interactions will be called ZBM-F and ZBM-0 interactions 
, respectively. The Thomas-Ehrmann shift is then taken into account through the combined 
effect of the mirror symmetry breaking continuum coupling and the modification of the s.p. 
energies of the effective SM interaction. 

The coupling matrix elements between the = 0^ g.s. wave function of ^^O and all 
considered states in ^^F and ^''O are calculated using the density dependent DDSMl inter- 
action |TO|JTT[] . The coupling to the continuum states is given by the matrix elements of 



'^75(L)/3 between g.s. of ^^O and all considered states in ^'^F and ^^O. The calculation of the 
radial wave functions and radial formfactors of the coupling to the continuum states goes 
as follows. The s.p. wave functions, which in the many body states J'^ of ^''F are not mod- 
ified by the selfconsistent correction to the finite initial potential, are calculated using the 
auxiliary potential f/('*"^). This potential, which contains central , spin-orbit and Coulomb 
parts, is adjusted to yield the binding energies of proton s.p. orbits 0(^5/2 and lsi/2 at the 
experimental binding energies of 5/2f and l/2f states in ^^F. This adjustement of binding 
energies of s.p. orbits makes sense because the many-body wave function 5/2^ (respectively 
l/2f) has a large amplitude of the component with one particle in 0^5/2 (respectively lsi/2) 
s.p. state outside of the ^^O core. The parameters of this potential are given in Table 1 for 
different values of the diffuseness parameter a. Without Coulomb term, the same potential 
is used also to calculate radial formfactors for neutrons in ^^F. In the calculation for ^^O, 
again the same potential with the Coulomb term of ^^O is used. Such a potential yields 
binding energies of neutron 0^5/2 and lsi/2 s.p. orbits very closely to the experimental bind- 
ing energies of 5/2f and l/2i states in ^^O. If the coupling between bound and scattering 
states modifies the s.p. wave function (pij in the many body state J'^, then the depth of 
the initial potential U{r) is readjusted to ensure that the energy of the s.p. state (pij in 
the selfconsistent potential U^'^^^r) is the same as the energy of this state in the auxiliary 
potential f/^'^"^). The remaining parameters : Rq, a, Vso, of the initial potential are the 
same as in iji'^'^^). This procedure for radial wave functions yields the same asymptotic 
property for a given s.p wave function in all different channels. 

The calculation of the first-forbidden /5— decays presented in this work are the extension 
of the calculations of Towner and Hardy [|T^ (see also Ref. |T^ for the presentation of 
the method). In the following we give only few elements of this approach to introduce 
the notation used in this work. The calculation of the absolute decay rate uses in,!^ : 
fti/2 = 6170 s, where ti/2 is the partial half-life and : 

f = C{W)F{Z,W){W^ -lY^^W{Wo-WfdW . (7) 

In the above expression, W is the /?— energy, Wq is the maximum /3— energy and Z is the 
charge of the final nucleus. The first-forbidden shape factor C{W) can be written to a good 
approximation as |I^p!2| : C{W) = k{l + aW + bW~^ +cW'^), where coefficients k, a, b and 
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c depend on the nuclear matrix elements, Wo and ^ = Ze^/2R. (For the nuclear radius R 

(8) 



we use the prescription of Wilkinson [15|.) Consequently : 

/ = k{Io + ah + + c/2) . 



The evaluation of integrals is given in [|T6|. Relation of the coefficients k,a,b,c to the 
nuclear matrix elements of different rank is |]T^,|T3 • 



(0) 



+ 



1 



1 



Ci + ^{x + uf- -fiiiMx + u) + —W^{2x + uf - Y^A2(2x - uf 



(1) 



+ 



+ 



12 



(2) 



ka 
kb 



1(0) 









.0 



l(2) 



(9) 



iiiii{-[Cowr + [Ci{x+u)f^} 



1 



kc = —[8^2 + (2x + uf + X2{2x - + 
18 



12 ^ 



(2) 



where 



Co = + ^wWo 



V = ^v + ^w , 

Y = ^'y-^{u' + x) , 

The quantities /xi,7i and A2 are defined in terms of the electron wave functions and have 



(i = Y + -{u-x)Wo 



values close to unity . The non-relativistic form of the nuclear matrix elements is [ff7 



w = AV3(J/r/|||ir[Ci X afT\\\JiTi)C 
X = -{JfTf\\\irCi ■ T\\\JiTi)C , 
u = XV2{JfTf\\\ir[Ci X a]^T\\\JiTi)C 
z = -2A(J/T/|||2r[Ci X a]'^T\\\JiTi)C 



(10) 



where A = -gA/gv = 1.2599(25) p[ and : 

(T,T,1±1|T^T,,) 



C 



[2(2J, + l)(2T; + l)]i/2 
The remaining matrix elements in the non-relativistic form are : 

^'v = -\V3{JfTf\\\^[cx X V]V|||J,T,)C 

C'y={JfTf\\\^v-T\Uh)c , 



(11) 



(12) 
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where M is the nucleon mass. It has been found by Warburton et al. p9|j20|] that the 
matrix element ^ f of the time-hke piece of the axial current is strongly enhanced by meson- 
exchange currents, mainly the one-pion exchange. The enhancement factor that multiplies 
the impulse- approximation axial-charge matrix element, has been determined by comparison 
to experiment for A ~ 16 to be e = 1.61 ±0.03. In all calculations discussed in this work, we 
multiply the matrix element ,^ f by a constant factor 1.61, like in the analysis of Borge et al. 
0]. Somewhat smaller enhancement factor has been used in the study of Millener [Q. This 
meson exchange contribution to the axial charge is very important, changing the absolute 
decay rates of the first-forbidden decays /~ and /"*" by a factor ~ 3.7 (see also the discussion 
in 0). However, its influence on the ratio f^/f~ is somewhat less important. 

The matrix elements w , x and u are obtained from w, x and u by including an extra 
factor in the radial integral []T^. For the higher energy part of the spectrum in light nuclei, 
the term kbW~^ can be neglected and the energy dependence of C{W) is determined entirely 
by the matrix elements of multipolarity 1 and 2 []T^, i.e. by ka and kc. As written above, 
the expressions for the matrix elements apply to electron emission. For positron emission 
we have to make the following replacement : Z —Z and A ^ — A. Even assuming perfect 
mirror symmetry, the nuclear matrix elements combine in / with different signs what gives 
a different shape correction factor for electron and positron decays. The evaluation of 
matrix elements in the basis of selfconsistently determined s.p. wave functions in the initial 
J'^ = 1/2^ and final JJ = 1/2+ many body states follows the procedure described by Towner 
and Hardy [O and adopted in Ref. [Rl. 



III. DISCUSSION 

An exceptionally large asymmetry of mirror /?— decays : ^^N(J'^ = l/2f) — >^'^0{J'^ = 
l/2f) and -'^''Ne(J'^ = 1/25") — >^'^F{J'' = l/2f), raises the question about a role of largely 
different radial sizes of lsi/2 s.p. wave functions in the initial state l/2j" of ^^Ne/"^'^N. In 
the final state l/2i , the large asymmetry of radial configurations of ^"^F, ^'^O and the halo 
structure of the ^''F is well known. A dominant component of the 1/2^^ wavefunction in 
-'^''Ne/-'^''N is 7r(0(i5/2)'^(0p7/2)/^(0'^5/2)^(0Pi/2) configuration outside the core of -'^^O but this 
configuration does not play a role in the first-forbidden /5— decay. The dominant contribution 
comes from the small component of the wave function 7r(ls^^2)^(0Pi/2)/^(-'-'^i/2)'^(0i'i/2) 
the g.s. of ^^Ne/^^N. For this component, lsi/2 nucleon makes a transition to fill the Opi/2 
hole. 

Assuming strict mirror symmetry in the /?— decays, the SM analysis yields the ratio 
/+//" which is much smaller (/'*'//" — 9.6) than the experimental value (24 ±4) 0. 
Allowing for the variation of radii of selected s.p. orbits, Borge et al. have shown that 
the observed values of the ratio /+//" can be reproduced in the SM analysis assuming 
a difference of ~ 0.6 fm between the oscillator length parameters for proton and neutron 
lsi/2 orbitals. It was assumed that radii of lsi/2 neutron orbits in ^"^N and ^^O are the 
same. Similar assumption has been made for radii of lsi/2 proton orbits in ^^Ne and in ^^F. 
Remaining s.p. orbits have been assumed to have the same radius in ^^F, ^'^Ne, ^^O and ^^N. 

High sensitivity of the /3— decay asymmetry to the spatial features of s.p. orbits [Q, 
provides a challenge for the SMEC approach because radial characteristics for certain or- 
bits as well as their asymptotic properties are determined consistently for each studied nu- 
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cleus and, moreover, model parameters ( the effective SM interaction, the residual coupling, 
radius/diffuseness of the initial/auxiliary average potential) are determined by analyzing 
different reaction and spectroscopic data in the same many-body framework. 

In Fig. 3 we show and /^/ /~, which are calculated in SMEC for different values of 

the diffuseness parameter of the initial potential. In each case, the same diffuseness is taken 
for initial potentials in ^^Ne, ^^N, ^^F and ^^O. The shaded areas show the experimental 
limits for these values. and /~ are calculated with ZBM-F and ZBM-0 interactions 
(the solid lines), respectively. The dotted line in the middle part of Fig. 3 shows results 
for /~ obtained with ZBM-F interaction. The corresponding ratio f~^/f~, shown by the 
dotted line in the lower part of Fig. 3, corresponds then to strictly mirror symmetric SM 
interaction. The coupling to continuum is given by DDSMl residual interaction in all studied 
cases. The calculations of the radial wave functions and radial formfactors of the coupling 
to the continuum states for ^''F/^^0 have been described above. In the calculations for 
^'^Ne/^'^N, we employed the initial/auxiliary potentials obtained from those for ^''F/^''0 by 
an appropriate modification of the Coulomb potential. One should mention, that initial 
potentials U{r) (t/(""^)(r)) for different diffuseness parameters, yield very similar spectra 
for each considered nucleus. Therefore, the energy spectra do not provide any constraint 
on the diffuseness of the average field. The experimental value of the ratio /"*"//" can 
be reproduced using U{r) ([/(»«^)(r)) with a very large value of the diffuseness parameter 
(a ~ 0.8 fm). However, none of these potentials with extremely thick skin can reproduce 
the experimental values for either /+ {f^xp = 927 ± 95) or /~ {fe^p = 44 ± 7). 

The overall contribution of matrix elements of rank 1 in /"*" and /~ increases with 
increasing diffuseness parameter from ~ 11% at a = 0.4 fm to ~ 21% at a = 0.8 fm. This 
strong relative increase is similar for and /~. More insight into the a— dependence of 
and /~ can be gained from Fig. 4 (see also Tables 2 and 3) which shows the matrix elements 
: klo (the upper part) and the sum kali + fcc/2, separately for ^''Ne(/3'^)^'^F (the solid lines) 
and ^^N(/5~)^^0 (the dotted lines) decays, klo is the combination of nuclear matrix elements 
of rank and 1. Coefficient k is the energy- independent term of the first-forbidden shape 
factor C{W). The sum of kali and /cc/2 depends on nuclear matrix elements of rank 1 
only. In our case, the nuclear matrix elements of rank 2 associated with the operator 
{c.f. equation (^) are identically equal zero. The coefficients ka and kc determine W— and 
W^^— dependence of the shape correction factor C{W). ZBM-F and ZBM-0 interactions 
are used in the calculation of j3~^ and j3~ decays respectively. {kIo)~^ and (klo)' decrease 
with increasing a and the decrease rate is similar in both cases. The strong increase of 
the ratio f^/f^ for large a is caused by a different a— dependence of {kali + kcl2)^ and 
{kail -\- kcl2)~ ■ With increasing diffuseness of U{r) (f/(°"^)(r)), {kali + kcl2)^ decreases 
strongly whereas {kali + kcl2)~ remains approximately constant and is negative. In the 
studied range of a values, the ratio {klo)^ / {kl())~ increases only by ~ 25% whereas [{klo)^ + 
{kail + kcl2)~^]/[{kIo)~ + {kali + kcl2)~] — f~^/f~ increases by ~ 70%. At the same time, 
kbl^i is small and almost constant. Hence, those nuclear matrix elements of multipolarity 
1 which determine the /3— energy dependence of the first-forbidden shape factor play also a 
salient role in exhibiting the asymmetry in the mirror first-forbidden /?— decays : ^'^Ne(/?+)^''F 
and ^^N(/3-)i70. 

This effect is independent of whether ZBM-F or ZBM-0 effective interaction is used 
in describing /3~ decay : ^''N(/3~)^''0. The dashed line in Fig. 4 shows the dependence of 
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{kIo)~ and {kali + kcl2)~ on the parameter a of the first guess potential U{r), which is 
calculated using the ZBM-F interaction, as in the calculation of ^'^Ne(/3^)^''F decay rate. As 
before, (fc/o)~ decreases strongly and {kali + kcl2)~ remains constant. This shows that the 
a— dependence of the ratio f^/f~ is not related either to the choice of the particular variant 
of ZBM interaction or to the assumption of strict mirror symmetry in (3^— and decays. 
On the other hand, the overall magnitude of this ratio depends strongly on the amount of 
mirror symmetry breaking in the many-body wave function. Replacing ZBM-0 by ZBM-F 
in the calculation of ^^N(/5~)^'^0 decay rate, leads to the decrease of /~ by a factor ~ 1.8 
(see also Fig. 4). This sensitivity, which has been stressed by Millener [Q, is stronger than 
the a— dependence of The separate contributions of (fca/i)^, are given 

in Table 2. Nuclear matrix elements are shown in Table 3 for different SM interactions and 
SMEC radial wave functions. /~ is calculated for both ZBM-0 and ZBM-F interactions. 
One can see that in the considered range of surface diffuseness parameters the enhancement 
rate of the ratio f'^/ f~ is ~ 1.4, independently of the assumption of strict mirror symmetry 
of the SM interaction. 

In order to quantify the relative importance of radial sizes of the s.p. states and the 
mirror symmetry breaking in the many body wave functions, one has to constrain the surface 
diffuseness of the initial potential. Such a constraint can be provided by the proton radiative 
capture cross section ^^0(p, 7)^'^F to the 'proton halo state' l/2f, which is very sensitive to 



the size of lsi/2 proton s.p. wave function ||11| , |21| . The astrophysical ^-factor as a function of 



the cm. energy is shown in Fig. 5 for three values of the diffuseness parameter of the initial 
potential : a = 0.4 fm (the dotted line), 0.55 fm (the solid line) and 0.8 fm (the dashed line), 
separately for each decay branch : ^'^0(p, 7)^'^F( J'^ = 1/2^) and ^^0(p, 7)^^F( J'^ = 5/2^^). 
One should notice that in all cases, energies of proton s.p. orbits 0^5/2 and lsi/2 in the 
selfconsistent potentials for diffenent a are at the experimental binding energies of 5/2f and 
l/2f states in ^''F, respectively. The sum of contributions from these two decay branches 
is shown in the lower part of Fig. 5. The calculations are performed using SMEC wave 
functions (5) and (6) which have been obtained exactly in the same way and for exactly the 
same input parameters as those used in the calculation of both the decay rate ^^Ne(/3^)^^F 
(see Figs. 3 and 4) and the spectrum of ^^F (Fig. 1). The scale of excitation energy is the 
same as cm. energy in the p O system. The photon energy is given by the difference 
of cm. energy of [^^O + p]jT system and the experimental energy of the final state (1/2^ 
or 5/2^) in ^"^F. We have taken into account all possible El, E2, and Ml transitions from 
incoming s, p, d, /, and g waves. It is clearly seen in Fig. 5 that the SMEC calculation 
for a = 0.4 fm underestimates the experimental capture cross-section. On the other hand, 
the calculation for a = 0.8 fm, for which the ratio f^/f~ agrees with the data (see Fig. 3), 
overestimates strongly the data. Realistic values of the surface diffuseness parameter, which 
are compatible with the proton capture data |2^, are : a = 0.55 ± 0.05 fm. In this range, 

agrees perfectly with the experimental data, whereas /~ is too big. For a = 0.55 fm, /~ 
overshoots the data by a factor ~ 1.4. 

Since the radial dependences which are consistent with the proton capture reaction data, 
give excellent fit of both the /J"*"— decay rate and the spectrum of ^''F, the discrepancy found 
for /~ and f^/f~ should be explained by the deficiency of the effective SM interaction 
to reproduce the configuration mixing in ^^O (^''N). If one includes the effect of mirror 
symmetry breaking through the modification of energies of s.p. orbitals of the SM, then 
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the experimental value of /~ can be reproduced by SMEC wave function for a = 0.55 
fm using : Sd^^^ = 3.21 {Sg-^/^ ~ 3-3). For this interaction, called ZBM-0*, Fig. 6 shows 
SMEC energies and widths for positive parity (l.h.s. of the plot) and negative parity (r.h.s. 
of the plot) states of ^^O. Negative parity states are reproduced better by ZBM-0* than 
by ZBM-0 SMEC calculations (c./. Fig. 2 and Fig. 6). The splitting of 5/2^ and 1/2^ 
states in SMEC/ZBM-0* is slightly larger than the experimental splitting. The spacing 
^<i5/2 ~ ZBM-0* is reduced too much as compared to ZBM-F interaction to account 

for the absence of the expected renormalization of the two-body matrix elements. Since the 
dominant configuration in 5/2f and l/2f states is the Ip-Oh component outside a closed 
core of ^^O with one nucleon in either or Si/2 shells, therefore it is quite natural that 
the splitting Sd^^^ — e^^^j ZBM interaction has a particularly strong effect on the relative 
position of these two states. 

An interesting supplementary quantity is the B{E2) transition matrix element between 
l/2i and 5/2^^ bound states of ^"^F and ^"^O. Assuming the effective charges : Cp = 1 + Sp, 
Cn = Sn, with the polarization charge 6 = 6p = 6n = 0.2, which are suggested by the 
theoretical estimates one finds in SMEC/ZBM-F with a = 0.55 fm : B{E2) = 74.85 
e^fm"^ for ^^F. The experimental value for this transition is B{E2)exp = 64.92 e^fm^. In ^^O 
one finds similar resuhs : B{E2) = 3 and 3.2 e^fm^ in SMEC/ZBM-0 and SMEC/ZBM- 
O* calculations, respectively. The experimental value for this transition is B{E2)exp = 6.2 
e^fm^. 

Fig. 7 shows /~ and /^/ /~, which are calculated in SMEC for different values of the 
diffuseness parameter of the initial potential. The shaded areas give the experimental limits. 

and /~ are calculated with ZBM-F and ZBM-0* interactions, respectively. An essential 
physical parameter is the amphtude of a component (ls^y20pr/2) S-^- wave functions 

of ^''Ne and ^''N. These amplitudes are shown in Table 4 together with the dominant Ip-Oh 
component of the wave function in the final state 1/2'^. One can notice that the amplitude 
of (ls2/20p-/2) configuration in ^^Ne (ZBM-F) is ~ 30% bigger than in ^^N (ZBM-0*). One 
should however recall that this important difference concerns the small component of the 
wave function and the large component of the many-body wave function in the g.s. of ^^Ne 
and ^^N differs by less than 5%. The same change of the effective interaction induces only a 
small (~ 12%) difference of the dominant Ip-Oh configuration in 1/2^^ wave function of ^''F 
and ^^O. 

In conclusion, we have found that the increase of the ratio f^/f~ is correlated with the 
increase of radius of weakly bound lsi/2 s.p. orbit in ^^F, in accordance with the conclusion 
of Borge et al. [0], but to obtain agreement with the data one has to assume an unrealistic 
geometry of the self-consistent potentials in ^"^Ne, ^''N, ^"^F and ^''O, which disagrees with 
the proton capture data. If, consistently, one takes into account the constraint from the 
capture data, one can estimate mirror symmetry breaking in the SM effective interaction 
in ZBM space for A = 17 nuclei. This effect is less than ~ 5% for the dominant term in 
the g.s. wave function and about ~ 30% for the small component of the wave function 
which plays a crucial role in the first-forbidden /?— decay. In the final nucleus, the mirror 
symmetry breaking is less than ~ 12%. These estimates are expected to depend somewhat 
on the SM space used. In particular, the absence of 0p3/2 and Oo?3/2 subshells in the ZBM 
space may lead to the amplification of the sensitivity in the lsi/2 Opi/2 contribution to 
the charge- dependent effects 0. 
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Many nuclear matrix elements contribute to the transition probability so the experi- 
mental determination of the lifetime and spectrum shape alone is usually insufficient to 
determine them all. A unique first-forbidden /5— transition from ^''N to the g.s. of ^^O is 
known (/ = 24 ± 8) [0]. Unfortunately, this transition tests the nuclear matrix element of 
rank 2 which are absent in the first-forbidden /5— transition and which cannot be calculated 
reliably in a small effective SM space. SMEC gives for this unique first-forbidden transi- 
tion a value which is ~ 3 times larger than the experimental value and depends weakly on 
the chosen hybrid of the ZBM force. This discrepancy is expected to be solved in calcu- 
lations using large basis up to 4huj [^,^. Only in favourable circumstances, the nuclear 
structure information can be unambiguously extracted from the first-forbidden /3— decays 
and relative importance of the configuration mixing (internal mixing) and the exotic radial 
dependences of s.p. wave functions (external mixing) can be disentangled. It seems that 
the first-forbidden /3— decays, which depend sensitively both on the fine details of the SM 
effective interaction and on the radial formfactors of the wave functions in mirror systems 
will be particularly difficult to exploit as a direct and unambiguous source of information 
about unstable nuclei. 
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TABLES 



TABLE I. The parameters of initial potentials U (r) used in the calculations of self-consistent 
potentials U'^^^\r) for si/2 and ^5/2 s.p. wave functions in 1/2^ weakly bound state and 5/2+ 
ground state of ^^F, respectively. The residual coupling of Q and P subspaces is given by the 
DDSMl interaction. For all considered cases the radius of the potential is Rq = 3.214 fm. For 
more details, see the description in the text. 
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0.50 


-53.6975 


1.5185 
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0.60 


-52.139 
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0.65 


-51.334 


5.475 


0.70 


-50.521 


6.706 


0.75 


-49.705 


7.898 


0.80 


-48.889 


9.054 


0.85 


-48.076 


10.178 
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TABLE II. Nuclear matrix elements and decay rates are shown for different values of the 
diffuscncss parameter of the 'first guess', initial potential U(r). Quantities for decays are 
calculated using ZBM-F interaction. For decay, ZBM-0, ZBM-0* and ZBM-F interactions 
are used. 
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TABLE III. Nuclear matrix elements for different values of the diffuscness parameter of the 
initial potential U{r). Quantities for decays are calculated using ZBM-F interaction. For 
/?- -decay. ZBM-0, ZBM-0* and ZBM-F interactions are used. 



Quantity 


Matrix element 


a = 0.4 


a = 0.5 


a = 0.6 


a = 0.7 


a = 0.8 


(ZBM-F) 




24.05 


22.62 


21.03 


19.43 


17.89 






-8.44 
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-6.78 


-6.21 
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-0.53 


-0.56 
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TABLE IV. The amplitudes of components {lsy2^Pi/^2^j ('-''^i/2^^'i/2) ™ initial state 1/2^ 
of the / P~ decay of ^^Ne/^^N and the dominant components lsi/2 and 0^5/2 in the final first 
excited state l/2f and in the g.s. of respectively, for different hybrids of the ZBM 

interaction. 



Interaction 
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0.41 
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0.69 
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0.399 


0.78 


0.665 


0.707 


ZBM-0 


0.311 


0.80 


0.62 
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ZBM-0* 


0.268 


0.81 


0.587 


0.637 
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FIG. 1. Comparison of experimental spectrum of F with the spectrum calculated using SMEC 
with the ZBM-F effective interaction. Residual coupling to the continuum state is provided by the 
density dependent DDSMl interaction [^] . The proton threshold energy is adjusted to reproduce 
position of the 1/2^ first excited state. The shaded regions represent the width of resonance states. 
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FIG. 2. The same as in Fig. 1 but for ^^O. The neutron threshold energy is adjusted to 
reproduce position of the first 3/2^^ excited state. 
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FIG. 3. /~ values and the ratio f^/f~ for the first-forbidden /?— transitions from the 

ground states of ^''Ne and ^^N, are calculated using SMEC for different values of the diffuse- 
ness a of the initial average potential. Calculations are performed in the ZBM space using the 
DDSMl residual interaction. The shaded areas show experimental uncertainty for /"*",/■" and 
/"'"//", respectively. The dotted line corresponds to the mirror symmetric SMEC calculations for 
^^Ne(/?+)^'^F and ^^N(/3-)i'^0 using in both cases the ZBM-F effective interaction. 
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FIG. 4. Nuclear matrix elements klo (upper part) and (kali + kcl2) (lower part) are plotted as 
a function of the diffuseness of the 'first guess' initial average potential. Calculations are performed 
in the ZBM space using the DDSMl residual interaction. The solid line shows the matrix elements 
for the transition ^'''Ne(/?+)^^F which is calculated using the ZBM-F interaction. The dashed and 
dotted lines show the matrix elements for the transition ^^N(/3~)^^0 calculated using ZBM-F and 
ZBM-0 interactions, respectively. 
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FIG. 5. The astrophysical S'-factor for the reactions ^^0(p, 7)^^F(J'^ = ^/^f) and 
^^0{p,'y)^'^F{J'^ = 1/2^), is plotted as a function of the center of mass energy Ecm for three 
different values of the diffuseness of the initial (auxiliary) potential : a = 0.4 fm (the dotted line), 
a = 0.55 fm (the solid line) and a = 0.8 fm (the dashed line). The experimental data are taken 
from Ref. [El. 
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FIG. 6. The same as in Fig. 2 but for ZBM-O* effective interaction which yields the correct 
value for the f3~ first-forbidden decay rate of ^^N in the ground state. The neutron threshold 
energy is adjusted to reproduce position of the first 3/2^ excited state. 
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FIG. 7. Z"*", / values and the ratio f~^/f for the first-forbidden transitions from the ground 
states of ^^Ne and ^^N, are calculated using SMEC for different initial values of the diffuseness 
of the initial average potential. ZBM-0* effective shell model interaction is used to calculate the 
structure of ^^O. For more details, see the description in the text and the caption of Fig. 3. 
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